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ABSTRACT 


A  simple  practical  numerical  method  for  estimating  the  wave 
resistance  and  the  wave  pattern  of  a  ship  operating  at  low  Froude  number 
is  presented.  The  method  is  based  upon  approximating  the  two-fold 
integral  over  the  ship-hull  surface  involved  in  the  usual  basic  expression 
for  the  far -field  wave-amplitude'  function  by  a  one-fold  integral  around 
the  ship  mean  waterline. 


ADMINISTRATIVE  INFORMATION 

This  study  was  funded  under  the  David  W.  Taylor  Naval  Ship  Research  and  Development 
Center’s  General  Hydrodynamics  Research  program.  Program  Element  61153N,  Project  Number 
SR02301.  Task  Area  SR0230101,  Work  Unit  1542-109. 
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INTRODUCTION 


A  numerical  method  for  evaluating  the  near-  and  far-field  wave  potential  and  the  wave  resistance 
of  a  ship  advancing  with  constant  speed  in  calm  water  within  the  theoretical  framework  of  the 
Neumann-Kelvin  theory  developed  in  Noblesse  [1]  was  recently  presented  in  Barnell  and  Noblesse  [2]. 
The  numerical  results  obtained  in  the  latter  study  show  that  a  very  large  number  of  panels  must  be  used 
for  approximating  the  ship-hull  surface  for  small  values  of  the  Froude  number.  The  numerical  method 
presented  in  [2]  therefore  becomes  less  practical  at  small  Froude  numbers,  as  is  indeed  noted  in  the 
conclusion  to  [2]. 

A  complementary  numerical  method  useful  for  small  values  of  the  Froude  number  is  presented  in 
this  study.  The  method  is  based  upon  approximating  the  two-fold  integral  over  the  ship-hull  surface 
involved  in  the  usual  basic  expression  for  the  wave-amplitude  function  by  a  one-fold  integral  around  the 
ship  mean  waterline. 

The  low-Froude-number  numerical  method  presented  in  this  study  can  be  further  simplified  and 
rendered  more  efficient  by  using  an  analytical  approximation  for  the  one-fold  integral  around  the  ship 
waterline  obtained  here.  This  analytical  approximation,  valid  for  sufficiently-small  values  of  the  Froude 
number,  is  presented  in  Part  2  of  this  study.  Numerical  applications  will  be  presented  in  Part  3. 


*A  complete  listing  of  references  is  given  on  page  17. 


BASIC  EXPRESSIONS  FOR  THE  FAR-FIELD 
WAVE-AMPLITUDE  FUNCTION 


This  study  considers  the  steady  potential  flow  due  to  a  ship  advancing  with  constant  speed  in 
calm  water  of  infinite  depth  and  lateral  extent.  Nondimensional  coordinates  x  =  X/L  and  flow 
variables  are  used,  with  the  length  L  and  the  speed  U  of  the  ship  and  the  density  c  of  water  selected  for 
reference.  The  mean  free  surface  is  taken  as  the  plane  z  =  0,  with  the  z  axis  pointing  upwards,  and  the 
X  axis  is  chosen  in  the  ship  centerplane  and  pointed  towards  the  bow. 

Equations  (42)  and  (32)  in  (1]  yield  the  following  expressions  for  the  nondimensional  wave 
resistance  r  =  R/^U^L^  and  the  nondimensional  velocity  potential  ^  =  <I>/UL  far  behind  the  ship: 

nr=  /  II  K(t)  ||2pdt,  (1) 

•'0 


n<j>{4)  ^ 


lm[E+(t:r)  +  E^ft;!)]  K(t)dt, 


where  E  +  (i;|)  represents  the  exponential  function 
E^ftif)  =  exp(v2p{p4  +  i(|±tn)}l; 
Furthermore,  p  is  defined  as 
P  =  (l+t2)''2, 


(2) 

(2a) 

(3) 


1/  is  the  inverse  of  the  Froude  number,  that  is 


V  !.r, 

and  K(t)  represents  the  far-field  wave-amplitude  function  defined  by  equation  (36)  in  [1]. 

I 

T  he  luneiion  K(t)  may  be  expressed  in  the  form 
Kit)  -  K  ^  (i)  +  K_(t), 


where  the  tiinelions  K  .  (t)  are  given  by 


K  .(t)  -  J^E  j(n^2  +  t^<|>,-n^ty<t>d  +  iv2p4,)  ty  dl 


4 


2/"  exp  (v2p2z)  E  .  (n^  f  v2p2^n  ^)da. 


(4) 


(5 


(6 


In  this  equation,  E  ^  represents  the  exponential  function 

E±  =  exp  [-iv2p(x±ty)J;  (7) 

Furthermore,  c  and  h  represent  the  positive  halves  of  the  mean  waterline  and  of  the  mean  wetted-hull 
surface,  respectively.  The  unit  vector  tangent  to  c  and  pointing  towards  the  bow  is  denoted  by  t(t^,ty,0), 
and  ir(nx,ny,n2)  is  the  unit  vector  normal  to  h  and  pointing  into  the  water,  as  is  indicated  in  figure  1. 

The  term  n  ^  is  defined  as 

n^  = -n2  +  i(nx±tny)/p.  (8) 

Also,  dl  and  da  represent  the  differential  elements  of  arc  length  and  of  area  of  c  and  h,  respectively. 
Finally,  represents  the  disturbance  potential  at  the  integration  point  iT  on  c  or  h,  represents  the 

derivative  of  4  in  the  direction  of  the  tangent  vector  t  to  c,  and  is  the  derivative  of  4  in  the  direction 
of  the  vector  nxt,  which  is  tangent  to  h  and  pointing  downwards  as  is  shown  in  figure  1. 

For  large  values  of  v,  or  more  generally  of  vp,  the  exponential  function  exp(v2p2z)  in  the  hull- 
surface  integral  in  equation  (6)  vanishes  rapidly  for  negative  values  of  z.  Therefore,  only  the  upper  part 
of  h  yields  a  significant  contribution  and  the  hull-surface  integral  can  be  approximated  by  a  line  integral 
along  c  in  the  low-Froude-number  limit.  This  low-Froude-number  asymptotic  approximation  of  the  far- 
field  wave-amplitude  function  in  terms  of  a  waterline  integral  is  obtained  in  the  following  section. 


Figure  1  -  Definition  Sketch  for  a  Single-Hull  Ship  with 
Port  and  Starboard  Symmetry 


LOW-FROUDE-NUMBER  APPROXIMATION  TO  THE  FAR-FIELD 
WAVE-AMPLITUDE  FUNCTION  IN  TERMS  OF  A  WATERLINE 

INTEGRAL 


The  mean  waterline  c  can  be  represented  by  the  parametric  equations 

x  =  Xo(A)andy  =  yo(A),  (S 

where  the  parameter  A  varies  between  its  bow  and  stern  values,  that  is  Ag^A^A^.  In  the  vicinity  of  c 
the  hull  surface  may  then  be  represented  by  the  parametric  equations 

x  =  X()(A)  +  x,(A)z  -I-  X2''A)z2  +  .  .  .  ( 

y  =  yo(^)  +  yiWz  +  +  •  •  •  ( 

where  and 

Equations  (7)  and  (lOa.b)  yield 

exp(v2p2z)E^  =  Eq±  exp(v2p2(i  ±)z] 

•  exp[-iv2p2(y2±z2-)-y3±z3+  ...)], 
where  Eq^  represents  the  exponential  function 
Eo±  =  exp(-iv2p2y^±), 
and  n^O,  is  defined  as 

yn±  =  (Xn±tyn)/p. 

The  Taylor  series  of  the  exponential  function  exp[-iv^%2^^^'*'y3*^^+  .  .  .  )]  is 

exp(-iv2p2(y2±z2  +  y^±z3+  )]  _ 

1  —  iv^p^z^T  * 

-v-p2z3[iy3±  +  v2p2z(y2±)2/2J  +  .  .  . 

The  parametric  representation  (10a,b)  of  the  hull  surface  yields 
da  =  in  dAdz, 

where  the  normal  vector  in  to  h  is  given  by 

in  =  (dx'dX)  X  id\/dz).  ( 

Equations  (lOa.b)  then  yield 

in  =  m„  +  m,z  +  m2Z^  +  .  .  .  ,  (16) 


where  the  vectors  m„  have  components 
tt^n  ~  ***11*) 

given  by 

n»n*  =  yn'*  *"11^=  -Xn'.  (16a.b) 

= ’‘o' yi-yo' *!•  (1^) 

=  Xj'  yi-yr  Xi  +  2(xo'y2-yo'x2).  (16d) 

>"2^  =  Xj'  yi-y2'  X,  +  2(x,'y2-y,'x2)  +  SCxo'yj-yo'xj).  (16e) 


In  equations  (16a-e),  and  hereafter  in  this  study,  the  superscript  '  denotes  differentiation  with  respect  to 
L  Thus,  we  have  x^'  =  dx^fAydA. 

The  potential  +(A,z)  on  h  in  the  vicinity  of  c  may  be  expressed  in  the  form 
MoW  +  +  4#)z2  +  .  .  .  .  (17) 

Equation  (IS)  yields 

(n^  +  v2p2^n  ^)da  =  (m^  +  v2p24m^)dAdz,  (18) 

where  we  have 

m^  =  -m^  +  i(mx±tmy)/p  (19) 

in  accordance  with  equation  (8).  By  using  equations  (16)  and  (17)  we  may  obtain 
m^  +  v2p2^rn^  =  v^pT^^mQ* 

+  (moX  +  v2p2z  +,^imo±)l 

+  z(m|’‘ +  v2p2z(,^^2± +<|'im|  ± +<|'2mo*)]  +  .  .  .  ,  (20) 

where  we  have 

m^i  =  -t- i(mn’‘±t  mny)/p.  (21) 

Equations  (20)  and  (14)  then  yield 

(m,(  + v2p2<^m  ^)  exp[-iv2p2(y2*z2  +  y3±z3  +  ...)) 

=  v^p^^Qmo*  +  Fi^  +  ^^2*  +  •  •  •  . 

where  we  have 

F,  ±  =  mo’‘  +  v2p2z(,^om,  ±  +<t'i  ^0"^) 

-iv‘»p-‘z2<}.()mo±y2*' 


(22a) 


Fji  =  m,*  +  v2p2z(^Qm2±++im,±+^2^0±-imo*y2±) 

-  iv4p'»z2rto(m ,  ±72*  +  >"0 *^3  * )  + +l"'*0^y2 *  1 

-v6p<iz  340^0=^  (72  *)^/2. 


Equations  (11),  (18)  and  (22)  yield 
exp(v2p2z)  E^(nj^  +  v2p2^n^)  da  = 

+f^l*  +zF2*  +  .  .  .  )  exp[v2p2(i  _iyj  ±)z]  dA  dz. 

We  then  have 


/exp(v2p2z)E^(n^  +  u2p2,|,n^)da  =  f  ^Eo±I±dA, 

A  Jxs 

where  I  ^  represents  the  integral 


r  = 


r 


(v^P^^qUIq  ±  +  Fj  ±  +  ZF2  ±  + 


)exp[v2p2(i  _iyj  ±)z]  dz. 


(22bj 


(23) 


(24) 


The  lower  limit  of  integration  in  the  integral  (24)  may  be  taken  equal  to  -  ».  We  have 
0 

z"exp[v2p2(i  -iy|±)z)dz  =  (-  l)nn'.  F2n  +  2q2n  +  2(u^)n  +  l, 


(25) 


where  q  and  u  ^  are  defined  as 

q  =  l/p  =  (l+t2)-l/2,  (26) 

u+  =  1/(1  -iy,±)  =  (1  -i(Xi±ty,)/pl->.  (27) 

Equations  (3)  and  (13)  were  used  in  equations  (26)  and  (27),  respectively.  By  using  equations  (22a,b)  and 
(25)  in  equation  (24),  we  may  obtain 

I  t  u  ±  +  F2q2u^G,  ±  +  F‘*q'Hu  i.)2G2±  +  .  .  .  ,  (28) 

where  G,  -  and  Gt*  are  defined  as 

-  '^y'o’  -  u±('t*omi=^+^imoM  -  2i(Uj^)2^omo±y2±.  (28a) 

G,"  =  -  y,'  +  2Ui(<t>om2* +'t'im,  ^  -'Vo'yz*) 

+  6i(U^)2[4iQ(m,±y2  *+"’()*  73*)  +  4'|'Tlo*y2*l 
-  12(u  .)3^omo±(y2M^- 


(28b) 


Equations  (17)  and  (10a, b)  yield 
4(X(,+x,z+  .  .  .  ,yo+yiz+  .  .  .  ,z)  = 

«l>(Xo-yo>0)  'I’jz  +  •  •  ■  • 


By  expanding  the  potential  in  a  Taylor  series  about  the  point  (XQ,yQ,0),  we  may  then  obtain 

•t*!  =  Xi4x  +  Vi'l'y  + -frz.  (29) 

where  <|>^,  represent  the  derivatives  of  <|i  with  respect  to  x,  y,  z.  We  have 

'I'x  =  V'j'-r  -  (i(>nn  +  (t>tr+<|ijnxT)*T; 

where 7  represents  the  unit  vector  along  the  x  axis,  and  4n  the  derivative  of  ^  in  the  outward  normal 
direction  n  to  h.  We  then  have 

K  =  ^  -  "z^y'^d-  (30a) 

We  may  similarls  obtain 

♦y  ^  'y*i  +  "z'x'^d*  (30b) 

'♦’/  -  »/'♦’!,  +  ("x'y  -  nytx)4’d-  (30c) 

l  et  D  be  defined  as 

1)  i,n^  t,n^  +  (X|t^  -  y,tjn^.  (31) 

Equations  (24),  (30a, b,c)  and  (31)  then  yield 
n  ('it,  *  y,t^)  -  <^1, 

where  the  relation  \|n^  +  y,ny  +  n^  =  n*ax(A,z)/az  =  0  was  used.  We  then  have 
M,'  *  I,*),  ^  n^2  +  (j>|n,ty/D 

‘  I's  (\,i,  .  ViiJn/^/DU,.  (32) 

On  the  mean  waterline  we  have  z  =  0;  equations  (15),  (16)  and  (16a, b,c)  then  yield 
n^  y,|  u(  1  t  £2) >  2^  (33a) 

n,  s„  11(1  f  £2)1/2,  (33b) 

11,  £(1^£2)12,  (33c) 

where  u  and  t  aie  detined  as 

“  K^n  )'  "  On' (34) 
f  ('•II  ''I  u.  (35) 

8 


Furthermore,  we  have 

dl  =  u  dA, 

‘x=-xo'/u.  ty=-yo'/u, 


(38a,b 


<t>,  =  --t-o'/u.  (35 

where  =  d4o(^)/dA.  Equations  (32),  (31),  (33a,b,c),  (34),  (35),  (36),  (37),  (38a,b)  and  (39)  then  yield 
+  =  -lyo'dA,  (4C 

where  I  is  given  by 

1  =  iv2p(^Q  +  Hi/u2(1  +£2)  (41 

with  H 1  defined  as 

Hj  =  (yo')2  +  eyo'u't'i 

+  [(1  +£2)xo'  -£yo'(x,Xo'  +y,yo')/uI.|.o'.  (42 

Equations  (40),  (7),  (9a, b),  (12)  and  (13)  then  yield 

h±  ("x^  +  di 

-'c 

fk 

=  -  /  Eo^Iyo'dA.  (43 

Equations  (6),  (23)  and  (43)  yield 


fh 

±(')-  J  Eo 


±A^  dA, 


where  A  ^  is  given  by 

A,  =  v2l  , -Ivo'. 

By  using  equations  (28)  and  (41)  in  equation  (45),  we  may  obtain 
Aj  =  v2(u^mo± -ipyo')(|>o-yo'H,/u2(H-£2) 

+  q2u^G,±  +  F2q4(u^)2G2±  +  .  .  .  . 


Equations  (27),  (21),  (16a,b,c),  (13)  and  (3)  yield 

u±mo±-ipyo'  =  =Fi(xo'±‘yo')(*/P^*yi)“±- 

Equations  (12),  (13)  and  (26)  then  yield 

v2Eo±(u^ino± -ipyo')'^o  =  ‘yi)“±+oJ' 

T  qEo±[(tqTiy,)u±^ol',  (47) 

where  the  superscript  '  denotes  differentiation  with  respect  to  the  parameter  A.  By  using  equations  (45), 
(46)  and  (47)  in  equation  (44)  we  may  obtain 


K  +  (t)  =  ±qlEo±(tqTiy,)u^+ol^S 


+  q- 


t- 


Eq  *a  ^dA, 


(48) 


where  a  ^  is  defined  as 

a±  =  u^ai^  +  F2q2(u j.)2G2*  +  .  .  .  ,  (49) 

with  a,  *  given  by 

a,±  =  +p[(tq7:iy,)Ui4ol'/Ui 

-  p2yQ'H|/u2(l +£2)u ^  +  G|±.  (50) 

At  the  bow,  we  have  y  =  0  for  z^O;  equation  (10b)  then  shows  that  we  have  yn(A)  =  0  for  A^Ag 
and  n^O.  Equations  (12),  (13),  (26)  and  (27)  then  yield 

±  qEo'Miq  +  iyilUi't'o  =  ±  tq2(exp( - iv2pxo)J<j>(/(i  “ 'Q^i) 
for  A=Ab.  Equation  (5)  then  shows  that  the  contribution  of  the  first  term  on  the  right  side  of  equation 
(48)  to  the  function  K(t)  is  null  for  A=Ab.  The  contribution  of  this  term  clearly  is  also  null  for  A=As  in 
the  usual  case  of  a  ship  with  a  closed  stern.  For  a  potential-flow-model  involving  an  open-stern  ship- 
and-wake  extending  to  downstream  infinity  we  have  as  and  the  first  term  on  the  right  side 

of  equation  (48)  therefore  also  vanishes  for  Equations  (5)  and  (48)  then  yield 

K(t)  =  q^y  ^(Eo  +  a^. +Eo-a_)dA.  (51) 
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By  using  equations  (42),  (28a),  (13),  (27),  (3S),  (21)  and  (16a-d)  into  equation  (SO)  we  may  obtain. 


after  some  algebraic  manipulations 

=  yo'A±/(l+£2)  +  2q(Xo' ±tyo')(u^)2B^4o 

+  C±<|>o'  +  uD^<|>,/(l  +£2)  +  ip(y ,<►()'  -yo'+i),  (52: 

where  A  ^ ,  B  ^ ,  C  and  D  ^  are  defined  as 

A±  =  [(l+pyo'/u)(l-pyo'/u)  +  £2]  +  i(pyo'/u)[y,(xo' ±tyo')/u  +  cl,  (53a; 

B±  =  q(y2=FtX2)  +  i(yiX2-x,y2),  (53b; 

C±  =  l(Tt-p2xo'yo'/u2)  +  (pyo'/u)2£(x,Xo'+y,yo')/u(l+£2)] 

+  i[yi(Xo'  ±tyo')/u+£][pxo'/u-(pyo'/u)£(x,Xo'+y,yo')/u(l+£2)],  (53c; 

=  ((Xo' ±tyo')/ul[(l +£2)(y,±iqt)Ui  +  i(pyo'/u)£y,J 

-  (pyo'/u)£(pyo'/u-i£).  (53d; 

In  the  particular  case  when  Xq'  ±tyQ'  =  0,  equations  (3),  (34)  and  (35)  yield  PXq'/u  =  -t, 
pyg'/u  =  ±  1,  (XiXg'  +yiyo')/u  =  “q(tXjTy,)  and  £  =  q(tyj  ±X|).  Equations  (52)  and  (53a,c,d)  then  yield 

a,±  =  ±  £[iyo' +(pyiTit)<t>o' Tu<|',]/(l±i£) 

+  ip(yi't’o' -yo'+i)  ifxg' ±tyo'  =  0.  (54] 

Furthermore,  if  the  hull  surface  intersects  the  free-surface  plane  orthogonally  at  the  point  (XQ,yo,0)  of 
the  waterline  for  which  Xg'  ±tyo'  =0,  we  have  n^^O.  Equation  (33c)  then  yields  £  =  0  and  equation  (54) 
becomes 

ai±  =  ip(yi4’o' -yo''t’i)  if  Xg' ±tyo' =0and  0^  =  0.  (55) 

Equations  (38a, b),  (39),  (29),  the  identity  <|'t  =tx*^x  +  *y*l’y  equation  (35)  yield 

(yi4o'  -yo'^l)/u=  + 

Equation  (55)  then  becomes 

a,  ±  =  iF2pyg'())^  if  Xg'  ±tyo'  =  Oand  n^  =  0.  (56) 

Equations  (49)  and  (56)  then  show  that  the  amplitude  function  a  in  equation  (51)  is  of  order  F2  at  a 
point  (Xg,yg,0)  where  the  phase  of  the  exponential  function  Eg±  =  exp[-iv2p(xg±tyg)]  is  stationary, 
that  is  where  Xg'  ±tyg'  =0,  if  the  hull  intersects  the  free-surface  plane  orthogonally  at  that  point. 


SUMMARY  OF  RESULTS 

The  mean  waterline  is  represented  by  the  parametric  equations 
x  =  Xo(A)andy  =  yo(A),  (5‘ 

where  the  parameter  A  varies  between  its  bow  and  stern  values,  that  is 

Ab  ^  A  Ag.  (57i 

In  the  vicinity  of  the  mean  free  surface,  the  hull  surface  is  represented  by  the  parametric  equations 

x=Xq(A)  +  X|(A)z  +  X2(A)z2  +  .  .  .  ,  (58i 

y  =yo(^)  +  yiWz  +  yzW^^  +  . . . ,  (58i 

where  Ag  ^  A  Ag  and  z  ^  0.  (58< 

The  velocity  potential  <t'(A,z)  on  the  hull  surface  in  the  vicinity  of  the  plane  z  =  0  likewise  is  expressed  in 
the  form 

+  .  .  .  .  (5! 

Differentiation  of  the  functions  x^fA),  yn(A),  i^nfA)  with  respect  to  the  parameter  A  is  denoted  by  the 
superscript  thus,  we  have  Xq'  =  dxgfAl/dA. 

For  sufficiently  small  values  of  the  Froude  number  F,  the  far-field  wave-amplitude  function  K(t) 
is  given  by  the  one-lold  integral 


(Eo  +  a^.  +Eo-a_)dA, 


where  we  have 


q  -  I  p  -  1  (I  +t2)l/2, 

and  I  „  •  and  a  .  arc  the  exponential  function  and  the  amplitude  function  defined  below.  The 
exponential  function  Eq-  is  given  by 
E()*  -  exp[  -  iv2p(Xo±tyo)), 
where  v  is  the  inverse  of  the  Froude  number,  that  is 


vvgi»  *■.•  yrv.K"' <- r-  <■ 

The  amplitude  function  a  ^  is  given  by 

1 1  1.^  ■mJt" W_"  «■- 

a±  =  u±ai±  +  F2q2(u^)2a2±  +  .  .  .  , 

where  u  ^  is  defined  as 

(64) 

u^  =.  1/[1  -iq(x,±tyj)], 

and  the  functions  aj  ±  and  a2^  are  now  defined. 

The  first-order  amplitude  function  aj  is  given  by 

=  yo'A  +  /(l-(-£2)  +  2q(xo' ±tyo')(u+)2B^^.O 

(65) 

+  C  +  uD^'^j/(l  +£2)  +  ip(yi'j»o'  -Vo'+i). 

where  we  have 

(66) 

u  =  [(xo')2  +  (yo')2j'/2, 

(67) 

t  =  (yo'xi-xo'yi)/u, 

and  the  coefficients  A  ^ ,  B  ^ ,  C  ^  and  D  ^  are  defined  as 

A±  =  [(l+pyo'/u)(l-pyo'/u)  +  £2] 

(68) 

+  i(pyo'/u)(y,(xo'  ±tyo')/u+£], 

(69a) 

B±  =  q(y2?tX2)  +  i(yiX2-x,y2), 

Ci  =  l(  +  t-p2xo'yo'/u2)  +  (pyo'/u)2£(x,Xo' +y,yo')/u(l +£2)1 

(69b) 

+  i  [y  ,(Xo '  ±  tyo '  )/u  +  £] [pxo '  /u  -  (py 0 '  /u)£(x ,xo '  +  y  ,y 0 '  )/u(  1  +  £2)] , 

-[(Xq'  ±tyo')/u][(l-i-£2)(y,±iqt)Ui+i(pyo'/u)£y,] 

(69c) 

-  (pyo '/u)«(pyo '/“-’«)• 

The  second-order  amplitude  function  a2*  is  given  by 

aji  =  -y,'  +  2u^(^o^2*+*i"’i^++2'"o*-'yo'y2*) 

+  6i(u^)2(^o(m,±y2*  +mo*y3*)  +  *l‘"o*y2*l 

(69d) 

-  12(u^)3^omo*(y2^)^- 

where  we  have 

(70) 

Xn*  =  q(Xn±tyn). 

(71) 

=  Mn  +  •q(yn'^‘’‘n')- 

(72) 

Ml  =  XiYi' -ViXj'  +  2(X2yo' -y2Xo').  (72b 

M2  =  xjyj'-yjXj'  +  2(x2yi'-y2X,')  +  acxjyo'-yjxo').  (72c 

At  a  point  where  the  phase  of  the  exponential  function  is  stationary,  that  is  at  a  point 


(XQ,yQ,0)  where  we  have  Xq'  ±tyo'  =0,  the  first-order  amplitude  function  ±  takes  the  form 

a,±  =  ±£[iyo'+(pyi±it)'i'o''pu'l‘il/(i±t£)  + 'P(yi4o'-yo''fri) 

if  Xq'  ±tyo'  =  0.  (73 

Furthermore,  if  the  hull  surface  intersects  the  plane  z  =  0  orthogonally  at  a  point  of  stationary  phase  we 
have 

a,±  =  iF2pyo'(t.^^  if  Xq' ±tyo'  =  Oand  n^  =  0,  (74 

and  the  amplitude  function  a^  then  is  of  order  F^. 

The  low-Froude-number  approximation  (60)  for  the  far-field  wave-amplitude  function  is  well 
suited  for  numerical  evaluation,  as  is  shown  in  the  next  section.  An  analytical  approximation  to  the 
integral  (60)  may  also  be  obtained  by  taking  advantage  of  the  rapid  oscillations  of  the  exponential 
function  Eg-  in  the  low-Froude-number  limit.  This  analytical  approximation  is  presented  in  part  2  of 
this  study. 

I  or  sut  ticiently  small  values  of  the  Froude  number,  it  may  be  acceptable  for  practical 
applications  to  approximate  the  velocity  potential  by  its  zero-Froude-number  limit,  in  which  the  free 
surface  becomes  a  rigid  flat  wall.  This  simple  approximation  for  the  potential,  together  with  the  low- 
Froude-number  approximation  for  the  far-field  wave -amplitude  function  obtained  in  this  study  and  the 
well-known  integrals  (1)  and  (2)  defining  the  wave  resistance  and  the  far-field  potential,  provide  a 
simple  practical  numerical  method  for  estimating  the  wave  resistance  and  the  wave  pattern  of  a  ship 
operating  at  low  Froude  numbers. 


NUMERICAL  EVALUATION  OF  THE  LOW-FROUDE-NUMBER  APPROXIMATION  TO  THE 

FAR-FIELD  WAVE-AMPLITUDE  FUNCTION 


The  functions  X|j(A)  and  yi((A),  where  k  =  0 _ 3,  in  equations  (58a, b)  can  be  determined  numerically 

by  considering  the  four  waterlines  z=  -kd,  where  k=0 _ 3  and  d  is  sufficiently  small  in  comparison 

with  the  (nondimensional)  ship  draft  d  (  =  D/L);  for  instance,  d  might  be  taken  equal  to  d/10.  Each  of 
these  four  waterlines  can  be  subdivided  into  an  equal  number,  say  n,  of  segments.  The  set  of  4(n-(- 1) 

nodal  points  (Xj,^,  yj^,  -  kd),  where  j  =  l_(n  +  1)  and  k  =0 _ 3,  define  n  -i- 1  constant  -  A  lines  A=Aj  in  the 

parametric  representation  (58a _ c)  of  the  upper  hull  surface.  For  a  given  value  Aj  of  A  corresponding  to 

any  one  of  these  constant  -  A  lines,  the  values  X|j(Aj)  and  y^fAj)  of  the  four  functions 

X|j(A)  and  yj^IA),  where  k  =0 _ 3,  can  be  determined  by  fitting  cubic  polynomials  through  the  points 

(Xj|(,  -  kd)  and  (yji^,  -  kd),  respectively,  for  k  =0_3.  Specifically,  the  cubic  approximation  to  the 
function  x(Aj,z)  takes  the  form 

x(Aj,z)~Xjo  +  (Xjo-Xj,)(z/d)  +  (Xjo-2xj,-t-Xj2)(z/d)(l -t-z/d)/2 
+  (Xjo-  3Xji  -I-  3Xj2-Xj3)(z/d)(l  -l- z/d)(2 -t- z/d)/6. 

A  similar  cubic  approximation  may  be  defined  for  the  function  y(Aj,z).  The  values  x,j(Aj)  of  the 
functions  X|^(A)  in  equation  (58a)  corresponding  to  the  foregoing  cubic  approximation  to  the  function 


x(Aj,z)  then  are  given  by 

Xo(Aj)  =  Xjo,  (75a) 

Xi(Aj)  —  ( 1  IXjQ  —  18xj ( -(- 9Xj2  —  2Xj3)/6d,  (75b) 

XilAj)  -  (2XjO-5Xj, +4Xj2-Xj3)/2d2,  (75c) 

X3(Aj)  -  (Xjo-3Xj, -(-3Xj2-Xj3)/6d3.  (75d) 


The  values  y|<(Aj)  of  the  functions  y|^(A)  in  equation  (58b)  are  defined  in  terms  of  the  y-coordinates  yjj^  of 
the  nodal  points  by  similar  expressions. 

The  values  '♦’k(Aj)  of  the  functions  ‘t'klA),  where  k  =0 _ 2,  in  equation  (59)  can  be  determined  in  a 

similar  manner  by  fitting  a  quadratic  polynomial  through  the  points  -  kd),  where  k  =  0 _ 2  and 

represents  the  value  of  the  potential  4  at  the  nodal  point  (Xj,^,  yji^,  -  kd).  This  quadratic  approximation 


is  given  by 


4(Aj,z)  4jo  +  (4jo“^jl)(z/‘5)  +  (4jo~2^jl +  +z/d)/2. 

The  corresponding  values  ^ij(Aj)  of  the  functions  4|j(A)  then  are  given  by 

=  'I’jo- 

=  (34j0-44ji  +  *j2)/2d. 

W  =  (4j0-2*j,  +  *j2)/2(52 


(76a) 

(76b) 

(76c) 


The  values  't'k(Aj),  where  k  =  0 _ 2,  and  the  values  X|j(Aj)  and  yk(Aj),  where  k  =0 _ 3,  determined  in 

the  manner  explained  above  for  a  set  of  n  +  1  values  Aj  of  the  parameter  A  can  be  used  as  base-values  for 
determining  the  functions  4k(^)»  ^k(^)  yk(^)  ®  denser  set  of  values  of  A,  for  instance  by  using 

cubic-spline  interpolation.  The  integral  (60)  defining  the  far-field  wave-amplitude  function  can  then  be 
evaluated  numerically  by  dividing  the  top  waterline  into  a  sufficiently-large  number,  say  N,  of  segments 
which  may  be  treated  as  straight  and  within  which  the  amplitude  function  a  ^  may  be  regarded  as 
constant.  Specifically,  the  integral  (60)  becomes 


K(t)~q2  2.  (Aj  +  ,-Aj)(aj+Ij+ +aj -Ij-), 
j  =  l 


(77) 


where  aj  -  represents  the  mean  value  of  the  amplitude-function  a  ^  within  the  segment  [Aj.Aj  ^  ,],  and 
Ij  *  is  the  integral  defined  as 


Ij  =  Ej  i  f  exp(  -  iv2p(6j  +  I  ±  t  )m)  d^, 
0 


with  Ej  ^  and  ±  defined  as 


Eji  =  exp( -iv2p0j±)and  0j±  =  X(,(Aj)  ±  ty,^A,). 

Likewise,  we  have 

Ej,,£  =  exp( -iv2p0j^,±)  and  0^,1 1  ,  ,)  t  ty„(Aj ,  ,). 

The  integral  Ij*  may  readily  be  evaluated  analytically,  with  the  result 

lj±  -  (Ej±+Ej  +  ,^)/2  if0ji  =  ^ 

=  iF2q(Ej,,±-E3±)/(0j^,i-0j*)  if0j^  0,  ,  ,  ^ . 


(78a,b) 


(79a,b) 


(80a) 

(80b) 


16 


The  number  of  segments  into  which  the  top  waterline  must  be  divided  increases  as  the  Froude 


number  decreases  [2].  The  foregoing  numerical  method  thus  becomes  less  efficient  as  F  becomes 
smaller.  However,  it  obviously  is  considerably  more  efficient  to  evaluate  the  (one-fold)  waterline  integral 
(60)  than  the  (two-fold)  surface  integral  involved  in  the  basic  expression  (6)  for  the  far -field  wave- 
amplitude  function.  Furthermore,  an  analytical  approximation  to  the  integral  (60)  may  be  obtained  by 
taking  advantage  of  the  rapid  oscillations  of  the  exponential  function  Eq*  at  small  values  of  the  Froude 
number,  as  was  already  noted.  This  complementary  asymptotic  approximation  is  presented  in  part  2  of 
this  study. 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 


1.  DTNSRDC  REPORTS.  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT 

2.  DEPARTMENTAL  REPORTS,  A  SEMIFORMAL  SERIES,  CONTAIN  INFORMATION  OF  A  PRELIM 
INARY,  TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION 

3.  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE  BY  CASE 
BASIS. 


